Effective Hamiltonian for transition-metal compounds. Application to Naa;Co02 
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We describe a simple scheme to construct a low-energy effective Hamiltonian Hcfs for highly 
correlated systems containing non-metals like O, P or As (O in what follows) and a transition-metal 
(M) as the active part in the electronic structure, eliminating the O degrees of freedom from a 
starting Hamiltonian that contains all M d orbitals and all non-metal p orbitals. We calculate 
all interaction terms between d electrons originating from Coulomb repulsion, as a function of 
three parameters {Fq, F2 and F4) and write them in a basis of orbitals appropriate for cubic, 
tetragonal, tetrahedral or hexagonal symmetry around M. The approach is based on solving exactly 
(numerically if necessary) a MO„ cluster containing the transition-metal atom and its n nearest 
O atoms (for example a CoOg cluster in the case of the cobaltates, or a CuOn cluster in the case 
of the cuprates, in which n depends on the number of apical O atoms), and mapping them into 
many-body states of the same symmetry containing d holes only. We illustrate the procedure for the 
case of Nai,Co02. The resulting //off, including a trigonal distortion D, has been studied recently 
and its electronic structure agrees well with angle- resolved photoemission spectra [A. Bourgeois, A. 
A. Aligia, and M. J. Rozenberg, Phys. Rev. Lett. 102, 066402 (2009)]. Affhough H^g contains only 
3d t2g holes, the highly correlated states that they represent contain an important amount not only 
of O 2p holes but also of 3d eg holes. When more holes are added, a significant redistribution of 
charge takes place. As a consequence of these facts, the resulting values of the effective interactions 
between t2g states are smaller than previously assumed, rendering more important the effect of D 
in obtaining only one sheet around the center of the Brillouin zone for the Fermi surface (without 
additional pockets). 



PACS numbers: 71.27.-fa, 74.25. Jb, 74.70.-b 



I. INTRODUCTION 

In condensed matter physics, several materials con- 
taining transition metal atoms (M) and oxygen have at- 
tracted great interest, for example the superconducting 
cuprates, manganites with high magnetoresistance, and 
more recently cobaltates and Fe pnictides. In the latter, 
the role of O in the electronic structure is replaced by 
either P or As. The electronic structure of these systems 
near the Fermi energy is determined essentially by some 
d orbitals of M and some 2p orbitals of O (or 3p, 4p 
for P, As in the pnictides), while the rest of the orbitals 
and elements play the role of reservoirs of charge. In 
most of these systems, strong correlations at M play an 
essential role. For example, while first-principles calcula- 
tions predict that the undoped cuprates like La2Cu04 
are non-magnetic metals^ they are in fact antiferro- 
magnetic insulators. In the cobaltates Naa;Co02, first- 
principles calculations^ predicted a Fermi surface with 
six prominent hole pockets along the T — K direction, 
which are absent in measured angle-resolved photoemis- 
sion spectrai^i^ 

Because of the presence of strong correlations, numer- 
ical methods which treat exactly the interactions in a 
small cluster of the system have led to a considerable ad- 
vance in the understanding of these materials.—"— Due 
to the exponential increase of the size of the Hilbert 



space, it is highly desirable that the set of relevant or- 
bitals included in these numerical studies is the smallest 
possible. For the cuprates, low-energy reduction proce- 
dures that eliminate the O degrees of freedom, simplify- 
ing the problem to an effective one-band one,— "— ^i^"— 
have been very successful, in spite of the fact that doped 
holes enter mainly at O atomsi^"— (optical properties re- 
lated with O atoms were calculated using these one-band 
models ^^■■^1-). Approximate treatments of the problem are 
also expected to lead to more accurate results when ap- 
plied to the effective model, because the construction of 
the latter begins with an exact treatment of the local in- 
teractions at AI. For example a slave-boson mean-field 
treatment applied to the three-band model (with Cu and 
O orbitals) for the cuprates predicts a metal for the most 
accepted parameters,— whereas the same treatment on 
the one-band effective model gives a gap that agrees with 
experiment 
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The main idea behind the construction of an effective 
one-band model for the cuprates, proposed by Zhang 
and RicefSS, starts from the exact solution of a CUO4 
cluster with one and two holes added to the vacuum in 
which the Cu atoms are in the Sd^*^ configuration and 
the O atoms in the 2p^ one. For two holes, the ground 
state is a singlet formed essentially by one hole in the 



3d.T, 



orbital of Cu and one hole in the linear combi- 



nation of 2p orbitals of the nearest-neighbor O sites, with 
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the same symmetry. Mapping this singlet state onto the 
corresponding state of a one-band model leads to the 
usual Hubbard or t — J models. Using a mapping in- 
volving non-orthogonal singlets, Zhang has shown that 
in a particular case, the mapping is exact. More sys- 
tematic derivations lead to additional terms of smaller 
magnitude,— i^i^ like three-site terms, which might be 
important for superconductivityi^i2Sri2^ Three-site terms 
play also an important role in effective models for vana- 
dates or other systems with t2g electrons.— Triplet states 
can be included perturbatively.— 

Second-order perturbation theory in the Cu-0 hopping 
applied to the three-band model for the cuprates leads to 
an interaction between Cu spins and O holes.— Adding 
to it the superexchange interaction between Cu spins one 
has the so called spin-fermion model. Similar derivations 
have been made for NiOr^ doped Y2BaNi05,— i2i and 
other perovskites.'^^ An important point is that the pa- 
rameters of the model are improved if instead of taking 
the values that result from perturbation theory, they are 
obtained fitting the energy levels of a CUO4 cluster in 
the case of the cuprates^ or a NiOe cluster in the case 
of the nickelatesi^ The Cu and O photoemission spectra 
for the cuprates obtained from the resulting spin-fermion 
model are practically identical to those of the three-band 
Hubbard model but required a much less computational 
effort for the same clusteri^ 

The above mentioned results suggest that in general 
the key for constructing a low-energy Hamiltonian is to 
solve exactly a cluster M0„ containing the transition 
metal atom and the O atoms of its neighborhood. In 
fact, the cell-perturbation method^" used extensively in 
the cupratesiSii^rJiii&iil exploits this idea: it divides the 
system into orthogonal cells containing one M atom each, 
using appropriate O Wannier functions centered at M, 
solves exactly each cell, writes the Hamiltonian in the 
basis of the resulting eigenstates, retains only the rele- 
vant low-energy eigenstates, and includes the rest per- 
turbatively. However, while in the cuprates the O Wan- 
nier functions can be easily constructed and the Hilbert 
space of the cluster involves only a few eigenstates, in 
the general case, as we show below, at least ten d and p 
orbitals are involved. In this case, the size of the Hilbert 
space of the M0„ cluster is 4^^^ and (except for favorable 
cases) the matrices should be solved numerically. For O 
rich systems, we propose to neglect the overlap between 
different O clusters and retain only the most important 
terms in the resulting effective model HcS (the effect of 
this overlap in the cuprates was calculated in Ref. [T5| ). 
A discussion of the effect of the overlap in the cobaltates 
is contained in Section [HI El Similarly to the above men- 
tioned ideas, the most relevant parameters are obtained 
from a fit of the low-energy levels of the original M—0 
multiband model H (obtained numerically) and those of 
HcS in the M0„ cluster. An additional calculation is 
required to obtain effective hoppings between clusters. 

An essential ingredient of the original multiband model 
H is the part of it which contains the interactions be- 



tween 3d electrons Hj . While it is straightforward to ob- 
tain it using known results of atomic physics, the deriva- 
tion is lengthy and it seems that only simplified forms 
were used so far for research in these systems. For exam- 
ple in a recent study of Fe pnictidesj^ a simplified ex- 
pression derived previously^ was used. While the form 
of Hi is well known when either only Cg orbitals'*'^ or 
only t2g orbitals^i^ are important, the correct expres- 
sions were not always usedj ^^'^^ In the general case, new 
terms appear which were not discussed before. One of 
the goals of this work is to present the complete Hj. We 
believe that this will be useful for future theoretical work 
on strongly correlated systems of transition metals, when 
the relevant orbitals cannot be restricted to either Cg or 
t2g only. 

In this work, we outline the derivation of HcS as de- 
scribed above. While for the sake of clarity we con- 
sider that the relevant non-metal orbitals are the 2p of 
O atoms, they could also be the 3p of P or 4p of As in 
the case of the Fe pnictides. We consider in detail the 
specific example of the cobaltates Naa;Co02. A brief de- 
scription of this case, together with a dynamical-mean 
field treatment of Hcs to study the electronic structure 
and Fermi surface of the system was published before4^ 

In Section II, we describe the construction of the multi- 
band model containing both M 3d (or 4d) and O 2p (or 
P 3p, As 4p) electrons, including all the d-d interaction 
terms. In Section HI, we derive the effective Hamiltonian 
HoS, which contains only effective d operators but no O 
ones. While some considerations are valid for the gen- 
eral case, the explicit construction is done for the case of 
Na2;Co02. Section IV is a summary and discussion. 



II. THE STARTING HAMILTONIAN 

As in many transition-metal compounds, we assume 
that the essential part of the electronic structure consists 
of the d electrons of the transition-metal atoms and the 
p electrons of the O (or other non-metal) atoms. An 
example in which only one d orbital is relevant is the 
three-band Hubbard model for the cuprateSf^i^ with pa- 
rameters determined by constrained-density-functional 
theory,—!^'* The Hamiltonian contains d-p hopping terms 
and interactions. The most important of the latter are 
the interaction terms among the d electrons. 

A. The interactions inside the d shell 

Here we can consider only one transition-metal atom 
and drop the site index for simplicity. The part of the 
Hamiltonian that contains the interaction among the 10 
d spin-orbitals is^ 

= ^ E ^AM-prf+d+d,d., (1) 
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where creates an electron or a hole at the spin-orbital 
A (Hj is invariant under an electron-hole transformation) 
and 



Vx^^p = J dridr2(px{ri)'fp{r2)- 



ri - r2 



(2) 

where ifxij^i) is the wave function of the spin-orbital A. 

In the basis of given angular momentum and spin pro- 
jections (A = 7n'^, cr'^), the orbital part of the wave func- 
tions can be written as R{r)Y™ (9 , (f)) , where R is the 
radial part and Yj^ is a normalized spherical harmonic. 
Using standard methods of atomic physics^i^ one ob- 
tains 



OO 

X J2 c'"{2m^, 2m'')c''{2m'', 2m'') • i^^ (3) 
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where P™(cos(0)) is a normalized Legendre function and 
r< (r>) is the smaller (larger) between ri and r2. The 
values of c*^ that are needed are tabulated in Ref. IH. 
To remove uncomfortable denominators, one defines the 
three free parameters as Fq = R'^, F2 = i?^/49 and F4 = 

In the presence of a cubic crystal field (point group Oh) 
or some other local symmetry (for example point groups 
Deh, Dih, Civ, Td), it is more convenient to change from 
the basis of operators with definite angular momentum 
projection dl^^ to that of irreducible representations of 
the point group (eg and t2g for Oh) orbitals using 



42. 



4i. = 



-'Oct 



d^ 



(5) 



In cubic symmetry, d^ 2 1 and dt o o correspond to 

J ji a: — y ,CT" 6z —r\a ^ 

eg symmetry and the remaining creation operators trans- 
form as the t2g irreducible representation. Each term of 
the resulting Hi contains two creation and two annihi- 
lation operators. Hi can be divided in four parts iJ„ 
according to the number n of eg operators present in 



each term. There is no term with only one t2g opera- 
tor and therefore H^ is absent. is usually enough 
to describe eg holes in late transition metals, such as 
nickelates,-- while Hq contains the relevant interaction 
terms for early transition metals with a few t2g electrons, 
such as titanates22. or rutenates<^ The sums over a (/3, 7) 
below run over the indices of eg (i2g) operators. The in- 
teraction can be written as 



Hi — H4 + Hq + Hi + H2 , 



(6) 



with 



H4 = 



+ (f/ - 2 Je) ^ 

(7i ,CJ2 

+ -^e ^ 42-y2 ,^j42-r2, 0-2 ^2;2 -2^2 0.2^32.2 _r2^o.j 

f 1 ,0'2 

+ -^e(42-j^2^^42-j^2,j_d3z2_r2,4.rf3z2_r2,t 



H.C.), 



(7) 



Ho - 



^ fj 2 J ^ ^ 



(8) 

Hi = XY [V3(42_,^2,,^d.x.CTi +H.C.) 

a I ,(T2 

X (4y,CT2'^a^,^2 +H.C.) 

- V3(42_j^2^^^rfy^,ai +^.C.){dly,^^d^^,^2 +H.C.) 

~ 2(422_^2^^^dxy,CTi +B..c.){dl^„^dy,^^„2 +H.C.) 

+ (42-r2,CTi^2^,'^l +H-C.)(4y,CT2^^a2.'^2 +H.C.) 

+ (42-r2,CTi^a^,^i +H.C.) 

X (4y,<.2^^-=^.-2 +H.C.)], (9) 

H2 = {U -2Jt) Y ^x'^-y^,aA^zx,a2 +»^yz,o-2) 

(T 1 ,CJ2 

^ ^ ^2;^ — ,0"! ^xy,cT2 

<7 1 ,0-2 

+ (t/-2Jb) ^ ^322_r2,CTi(^2a;,o-2 + ^^-^2,^72) 

CT"! ,C2 

+ (C/-2Je) E 

'^3z2 -,.2 n'xy,<T2 

Cl ,<T2 
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(71 ,(72 

^ '^l:^-y^,(7i'^yz,(72'^x^-V^,(72dyz,(7i) 

cri,cr2 

+ ^^6 '^L2-r2,(Ti'^=cy,o-2'^32:^-r2,CT2'^2;a,CTi 

17 1 ,0-2 

+ '/t['^^2_j^2 .|.C?^2_j^2 J^C'^ZX.i'^ZX,! + dyz,ldyz,^) + H.C.] 

+ >^a[c?^2_j^2 -[.'^^2_j,2 j_rfa:j/,4.'^a:j/,t + H.C.] 

+Jb[dl^2 ^^d^3^2 _r2 ^^{dzx,idzx,f + dyz^idyz^f) + H.C.] 

+ Je ^ ['^L2-r2,t'^L2-r2,4,^a:y4da:j/,t + H.C.] 

fl ,0"2 

+ 2A ^ (n-yz,<Ti - "-ZX,<Ti)(422_^2^0.2dx2-a^O-2 +H.C.) 

cri,a-2 

+ (4z2-r2,(Ti'^z2;,<T2'^2:^-J/^.cr2'^z2;,(Ti 

(71-72 

~ 4z2-r2.<Ti4z>'^2'^2;2-j,2 (j^C^J^Z.O-i + H.C.) 

+A[(4,t4,i-4z,t<z,i) 

X (da;2_2^2 J_d3^2_^2^^ - d2,2_j^2_^d322_^2 J^) + H.C.], (10) 

where U = F0+AF2 + 36F4, Je = AF2 + I5F4, Jt = 3^2 + 
2OF4, Ja = 35F4, Jt = F2 +3OF4, and A = \/3(i^2 - SEi)- 
The largest energy in Hi is the intraorbital repulsion 
U. Most terms of Hj (in particular all those of the pure 
Cg part Hi and the pure t2g part _ffo) involve two and only 
two orbital indices, like the (four different) ferromagnetic 
exchange interactions Ji, of the form 

E 4,^A.^2di,72d^,7i = -(2Se • S, + 1/2) (11) 

(71 ,72 

where — J2xx' 4 x'^^^"^i'^' ^'^^ ^P^^^ orbital 
^. The interorbital repulsion between these orbitals is 
related to the former by [/' = [/— 2 due to the spherical 
symmetry of Hj. For the same reason, = Ji, where 
Jj' is the energy for transfer of intraorbital pairs (like 
dp .^d'p ^dj^^dj^-f). The remaining terms, with prefactor 
proportional to A, involve more than two orbitals. 

The parameters F2 and F4 related to the Ji and A are 
expected to be very weakly screened in the solid as com- 
pared to the free atom and have little variation among 
the 3d series. For example, a fit of the lowest atomic 
energy levels (given in Rcf ''^) of V^"^ and Ni give nearly 
the same values F2 — 1300 cm~^ — 0.16 eV, and F4 — 88 
cm~^ = 0.011 eV within 2%. This results in an exchange 
energy Jg = 0.81 eV for Cg electrons and Je = 0.70 eV 
for t2g ones. This implies for example that for two Cg 
holes in cubic symmetry (as in Ni"*"^), the triplet ground 



state is separated by the excited singlet by 2Je ~ 1.6 
eV if covalency can be neglected. Recent calculations in 
3d metals suggest that the exchange interactions are re- 
duced in 30% in comparison with the atomic values4^ 
This reduction is also assumed in Ce compounds4^ 

In contrast to F2 and F4, Fq which determines the 
intra-orbital repulsion U is significantly screened in the 
solids and is difficult to determine theoretically. For the 
cuprates ?7 ~ 10 eV has been estimated by constrained- 
density-functional calculations^!^ and it decreases to 
the left of the periodic table inside the 3d series, because 
the 3d orbitals are more extended as a consequence of the 
smaller nuclear charge. The value of U can be extracted 
from optical experiments. 



B. The full M-O Hamiltonian 

Since usually there are only a few O holes present, it 
turns out to be more convenient to write the Hamiltonian 
in terms of hole operators (which annihilate electrons) 
acting on the vacuum state in which all transition-metal 
(M) ions are in the d^° configuration and the O (or P, 
As) ions are in the p^ one. The most important of the 
remaining terms of the starting Hamiltonian are the Co- 
O {ig^ below) and 0-0 hopping (t^^), parameterized as 
usual, in terms of the Slater-Koster parametersi^. We 
include a cubic crystal field splitting etg^ — egg = WDq 
at the metal sites. Extension to tetragonal, tetrahedral 
or hexagonal crystal fields is straightforward, while other 
symmetries may require a change in the chosen basis for 
the d orbitals. 

The Hamiltonian takes the form 

H= E '^e.dl^^d.aa + E '^*2g 4/3<Tf^»/3<^ 
i.a^e-g.a i,P&t2g,(7 

+ J2^^Pjv<yPr>T^+ E i5Hpl+s,nAi(7+ii-C-) 

+ E r'k>Lp^^^ + T.H} (12) 

Here p^^^. creates a hole on the O 2p orbital 77 at site j 

with spin a. The operator dl^^ has an analogous meaning 
for the M d orbitals at site i. The interactions at this 
site Hj has the form of Eqs. (jH]) to ((TU)l with the site 
index i added to the subscripts of the d operators. The 
subscript i + 6 in pl_^_g labels the different O atoms in 
the immediate neighborhood of the M atom at i (their 
nearest neighbors in highly symmetric structures). 

We have neglected here for simplicity the on-site 
Coulomb repulsion at O sites Up and the M-0 inter- 
atomic repulsion Upd- Experience in the cuprates indi- 
cates that the former is not very important, while it 
complicates the numerical treatment of the basic M0„ 
cluster (see below). Upd can be incorporated easily and is 
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important in the formation of excitons^i^ and in charge- 
transfer instabihties^^ which are beyond the scope of this 
work. 



III. CONSTRUCTION OF THE EFFECTIVE 
HAMILTONIAN 

A. Diagonalization of one cell 

The solution of the basic Af 0„ cluster containing a 
transition-metal atom and the O atoms in its neighbor- 
hood is greatly simplified for Up — 0, since among the 
3n O orbitals, only some linear (bonding) combinations 
with the same symmetry as the d orbitals hybridize with 
them, while the remaining (non-bonding) orbitals decou- 
ple. The presence of Up would introduce scattering be- 
tween the bonding orbitals and the non-bonding ones and 
does not modify the essential physics,— particularly for 
low or moderate number of O holes. 

For a perfect MOg octahedra (point group Oh) with M 
at the origin of cooordinates and the O atoms at 5 = ±x, 
±y, ±z, there are 13 O 2p non-bonding orbitals and the 
5 bonding O 2p orbitals are (spin indices and site index 
i are omitted for simplicity) 





2 iP'^jX P—yi,x 


Py,y 


+ P- 


-y.y)' 


PSz^-r^ = 




-z.z 








Px,x ~l~ P—x,x ~ 


Py^y ' 




y,a)' 


Pxy — 


1 


" Py.x 


-P 


-y,x)j 


Pyz = 


1 , 


" Pz,y 


-P- 




Pzx = 


1 

2 \Pz^x P—z,x ^ 


-P^,z 


-P- 


-X,.) (13) 



Writing the hopping terms in this basis, and using the 
Slater-Koster formulas^ the bonding part of the Hamil- 
tonian ([T^ in the MOq cluster takes the form 



-V3(pda)(dLpa^ + H.c.)} 

+ i'^t2,dl^dp„ + i^O + '2tp]plaPpa 

+2{pdTT){dl^Pf3^ + H.c.)} + Hi, 



where in terms of Slater-Koster parameters tp = 
-[{ppa) - {pp7t)]/2. 

If the cluster is elongated along the z direction or O 
atoms at ±z are missing, in general one should consider 
two different 2p orbitals that hybridize with (i322_,.2 . The 
first can be chosen as {—px,x + P-x.x ~Py,y +P-y.t/)/2, 



while the second is {pz.z — P-z.,z)l\/^ for D^h symmetry 
with obvious changes in absence of one or both O atoms 
at ±z. Similarly two independent states that hybridize 
with dyz in D^h symmetry are {py^z — P-y.z)/\/2 and 
(Pz.jf ^ P~-z.,y)l and the same changing y by x. The 
hopping terms that involve these (at most 8) O orbitals 
can be easily constructed using Ref. [sO- 

In any case, the size of the Hilbert space of iJb with 
at most 13 X 2 spin-orbitals is small enough to allow us 
to obtain its low-energy eigenstates numerically by the 
Lanczos methodi^ 



B. The cobaltates 

For the case of Naa;Co02, we have solved numeri- 
cally a CoOg cluster using the Hamiltonian given by 
Eq. mapped the corresponding states into those 

of an isolated Co atom with the corresponding charge 
to calculate effective on-site interactions, and calculated 
the hopping between effective Co sites mediated by O. 
The parameters of i/h were taken from fits to optical 
experiments^ as described below. We have neglected 
the trigonal distortion and assumed Oh symmetry in the 
cluster. As shown above, this assumption reduces the 
size of the relevant Hilbert space and in addition simpli- 
fies significantly the mapping between states of Ji and 
ffoff, because as we shall see, the symmetry identifies un- 
ambiguously the correspondence between states of both 
Hamiltonians. The effective trigonal crystal-field split- 
ting A = 313 '-^ 0.3 eV was given by quantum-chemistry 
configuration-interaction calculations,^^ As we discuss in 
more detail below, this is one order of magnitude smaller 
than the effective cubic crystal-field splitting and does 
not affect our main findings. 

Since Co^"*" is in a 3d^ configuration (3 holes in the d 
shell), the states of the CoOg cluster with n -f 1 holes 
are represented by a Co""*" ion in 7?eff- The most rel- 
evant values of n for Naa;Co02, are 4 (the formal va- 
lence of Co for a; = 0) and 3 (formal valence for the fully 
doped compound), but it is important to consider also 
n = 5 to calculate the effective interactions, as we shall 
show. As the vacuum at one site for the effective Hamil- 
tonian with only Co sites, it is convenient to chose the 
Co'^+ 3d^ configuration, occupied with the four holes. 
Thus Co''+ has one t2g hole and Co^+ has two t^g holes. 
Therefore, one expects that the interacting part of _ffoff 
at each site has the same form as Hq [Eq.®] but now 
t/^ff ^ J/cff ~ 2 Jeff , and J^g ^ Jeg because the cubic crys- 



tal field 



lODg reduces the symmetry from 



(14) that of the full rotational group to Oh'- 



Hf = U^^Y. 



,a<7' 
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TABLE I: Eigenstates and energies of one site of -ffcft , corresponding energies in the CoOe cluster and distribution of holes in 
the latter. For degenerate representations only one state is shown. 



Symmetry 


eigenstate 


E - 




E (eV) 


d Eg 


d t2g 


P Sg 


P t2g 











8.141 


3.03 


0.02 


0.94 





Tig 









17.271 


2.67 


0.91 


1.29 


0.14 


Alg 

E° 

Tl 




uu 


— JcS 
+ JcS 
~ JcB 


28.609 
28.096 
28.028 
27.320 


2.41 
2.40 
2.39 
2.34 


1.46 
1.51 
1.54 
1.63 


1.48 
1.57 
1.56 
1.65 


0.66 
0.52 
0.51 
0.37 



(15) 



eo = 13, et2„ = 1-2, te„ = 0- 



(16) 



This is the efFcctive Hamiltonian at one site, except for 
an unimportant eonstant C and chemieal potential term 
-/iiV, with N = S/3cr'^/3'^- tilde above the oper- 

ators is to remind us that the effective operators 
entering Eq. (jTS]) are different from those entering the 
starting Hamiltonian iJ, as discussed in more detail be- 
low. The resulting eigenstates and energies of Hcs at 
one site, are listed in the first three columns of Table I. 
It turns out that the low-energy eigenstates of the CoOe 
cluster with 4, 5 and 6 holes are well represented by the 
corresponding eigenstates of Hf^ . The fourth column 
of Table I corresponds to the lowest energy levels of the 
cluster Hamiltonian iJt, in each symmetry sector, using 
parameters determined previously by us^"* from a fit of 
polarized x-ray absorption spectra of Na2;Co02 to the re- 
sults of a CoOe cluster with 4 and 5 and holes including 
the core hole. In comparison with the present calcula- 
tions, the previous ones were simplified neglecting the 
exchange, transfer of intraorbital pairs, and terms pro- 
portional to A of Hi and H2, because their magnitude 
is smaller than the energy of the crystal-field excitations 
related with these terms. The effect of hybridization in- 
creases the splitting between t2g and Eg orbitals to more 
than 3 eV ^ F2, Fir-'^ The effect of the neglected terms 
is discussed in the next subsection. The parameters of 
iJfa in eV are 



The choice of the origin of on-site energies e^g = is 
arbitrary. The resulting value oi U — 4.516 eV and the 
charge transfer— energy Act — 2.9 eV are similar to 
those derived from other x-ray absorption experimentsi^ 

The most relevant eigenstates of the cluster for the 
physics of the cobaltates are those corresponding to the 
first two rows of Table I, corresponding to formal Co'^+ 
and Co^+ respectively. These states are separated by 
at least 2 eV from other eigenstates that cannot be rep- 
resented by i?off, like those containing non-bonding O 
orbitals or states with intermediate spin.™ This energy 
difference is one order of magnitude larger than the cor- 
responding hopping amplitude, assuring the validity of 
H'^^ as representative of the low-energy physics of H. 
We note that while in Co compounds, intermediate or 
high spin states are usual, experiments show that the 
low spin states of Co are present in Naa;Co02r^'^""— in 
agreement with our results. This is due to the fact men- 
tioned above, that in this system the splitting between 
t2g and eg orbitals (due to IQDq and Co-0 hopping) is 
more than 3 eV, considerably larger than the exchange 
energies. 

From the six energies listed in the fourth column of Ta- 
ble I and their corresponding expression for Hf^ listed in 
the third column, we have determined the four parame- 
ters of Eq. ([T5|) and the irrelevant constant C and shift 
in chemical potential /is- The result in eV is 



Fo = 3.5, F2 = 0.2, F4 = 0.006, 

— \/3 

(pdi:) = — - — {pda) ^ l,tp — 0.5, 



C/cff = 1.865, = 1.272, Jeff = 0.354, J^^ = 0.171 

(17) 



These parameters are smaller than those assumed in different calculations of the electronic structure of 
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Naa;Co02 near the Fermi surface which include the 
effect of correlations. ^^-S, For example values of Ucs 
between 3 eV and +00 were used. For our moder- 
ate values of the correlation energies, the hole pock- 
ets are still predicted in theory in contrast to experi- 
ment if the small first-principle value of the trigonal dis- 
tortion D is assumed.'fiiiS Instead, using the value of 
D given by quantum-chemistry configuration-interaction 
calculations these pockets are absent and the elec- 
tronic dispersion near the Fermi energy agrees with 
experiment!^ 

The reduction in the effective value of UeS compared 
to the 3d value U ~ 4.5 eV is due to screening effects of 
the full model that are "hidden" in iJoff. Actually, the 
operators d^^ of Hcs do not correspond to pure 3d holes, 
but are complicated creation operators which involve 2p 
states. The detailed expression of these effective oper- 
ator in terms of those of H is beyond the scope of the 
present work. Examples of the construction of effective 
operators are given in RefsJ^i^ii^. In this context, it is 
interesting to note that charge and current operators are 
expressed as pure spin operators when Heisenberg-like 
models are used as effective Hamiltonians for Hubbard- 
like models.— In any case, it is easy to understand 
that if the effective 3d operators d^^ have an important 
component of O 2p orbitals distributed in the cluster, 
the effective repulsion between electrons occupying these 
states is smaller than that between the corresponding op- 
erators localized on the same Co ion. 

In addition to covalency, another important point re- 
lated with it and noticed before by Marianetti et al.J'^ 
is that the addition of a new effective hole causes a re- 
distribution of the remaining charge. The distribution 
of charge among the different orbitals in the cluster for 
the different states is listed in the last four columns of 
Table I. The vacuum |0) of Hcs at one site, represents 
the CoOe cluster with 4 holes. Near 3 of them occupy 
the 3d Cg states and the remaining one occupies mainly 
linear combinations of 2p O orbitals with Cg symmetry 
[corresponding to the first two lines of Eq. ([T^]. When 
an effective 3d t2g hole is added, nearly 0.35 3d eg holes 
are promoted to the O states with the same symmetry. A 
similar effect is caused by the addition of the second ef- 
fective 3d t2g hole (with larger magnitude for the states 
of lowest energy). This is a combined effect of the in- 
terorbital Coulomb repulsion, which increases the energy 
of the Cg holes when t2g holes are added, and the strong 
hopping between 3d Cg holes and linear combinations of 
O holes with the same symmetry. 

The degree of covalency of the t2g states is less than 
that of the eg ones, due to the fact that the 3d-2p hopping 
is a factor 2 larger for the latter [see Eqs. (|14l) and p^ ]. 



C. The effect of Hi and H2 

Since most of the terms involved in Eqs. ^ and PU)) 
are cumbersome, one might ask what happens if one ne- 



glects them. The first four terms of H2 (the simplest 
ones), containing the Coulomb repulsion between d eg 
and t2g electrons are crucial to obtain the right distri- 
bution of particles and should be retained. Keeping the 
prefactors and neglecting the rest of the terms of Hi and 
H2, the cubic Oh symmetry is lost. Therefore, to re- 
tain the Oh symmetry, we have replaced these Coulomb 
repulsions by the average of all of them. 

All energies increase with respect to those of the full 
Hamiltonian. In particular for the states corresponding 
to nominal Co^"*" (second row of Table I), the energy in- 
creases by 0.92 eV. However, the distribution of holes is 
not dramatically affected. There is a 10 % reduction in 
the amount of d t2g from 0.91 to 0.82 which distribute 
evenly among other orbitals. Similar effects take place 
for nominal Co^^. The resulting effective repulsions are 
decreased considerably. Proceeding as indicated above 
we obtain Uctf = 1.47 eV, C/^g = 0.84 eV (35% less than 
for the full Hamiltonian), Jeff — 0.39 eV, and J^g — 0.25 
eV. 

Of course, the changes should be more dramatic when 
the effective cubic splitting {lODq plus effects of Co-0 
hopping) is smaller (like for other Co compounds), since 
for spherical symmetry all terms of exchange and transfer 
of intraorbital pairs are equally important. 



D. The effective hopping 

As shown first by Koshibae and Maekawaf^i the main 
hopping path from one Co site to a neighboring Co 
site in Naa;Co02 is via an intermediate O site. See 
Fig. [1] The expression of this effective hopping in 
second-order perturbation theory in the Co-0 hopping 
is toff = (pc?7r)^/AcT- Using the parameters given in Eq. 
p6)) . tcs = 0.34 eV is obtained. Clearly, this value is an 
overestimation, expected in covalent systems for which 
the hopping term is not small enough compared to the 
charge-transfer energy. An estimation based on the band 
width obtained from first principles gives icff = 0.1 eVM 

It is known that the cell-perturbation method allows 
us to obtain more accurate values of the effective hopping 
involving a calculation of only first order in the pertur- 
bation, using the eigenstates of the cellJ^iiSi^iiii^iil Here 
we explain this procedure for Naa;Co02, using the eigen- 
states of the CoOg cluster with 5 holes (represented by 
one occupied t2g orbital in i/cff) and 4 holes (the vacuum 
in iJoff). 

To be specific, we start from a state in which the cluster 
i (at the left in Fig. [T|) is in its ground state \gi{5,yz)) 
with 5 holes and yz symmetry and a neighboring clus- 
ter j (whose center corresponds to the rightmost orbital 
represented in Fig. [1]) is in its non-degenerate ground 
state 15^(4)) with 4 holes (symmetry Aig). After the 
hopping the final state is \gi{4)gj{5, zx)). The initial and 
final states are not orthogonal, but as we show below, 
this does not affect significantly the resulting effective 
hopping. The spin is conserved and the corresponding 
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FIG. 1: (Color online) Scheme of the main orbitals involved 
in an effective hopping from a dyz orbital at the left to a dzx 
orbital at the right through an intermediate p^ orbital. 

subscript is dropped for simplicity. This process is rep- 
resented in Hcs by the hopping dly^ ,^ — dl^^^.a- By 
symmetry, only the hopping terms of Eq. (jl2p which 
annihilate a hole in the 2pz orbital of the common O 
atom {pi+y^z = Pj-x,z, middle orbital represented in 
Fig. [T]) contribute to teg. Using Eq. with Slatcr- 

Koster parameters, Ea. (IT3l) and symmetry one has (call- 
ing tpd = (pdir), tp = -[{ppcr) - (pp7r)]/2) 

- toff = {g^{i)gj{5, zx)\{tpdd] ,^^ + tp{p]^^^^ - 
xpj-x,^|.9»(5,?;2;)5j(4)) 
= {9j (5, zx) I - tpdd]^^^ - tppl^^ \gj (4)) 

x(.9»(4)|ip,.|g.(5,yz)) (18) 

In the equation above, the part of the Hamiltonian that 
leaves the initial state \gi{5,yz)gj{4:)) unchanged is sub- 
stracted [see Eq. ([22l) below]. We have checked numeri- 
cally that the perturbative expression tcs = (pdir)'^ / Act 
is recovered for large Act- 

To calculate the matrix elements, we have used the 
simplified Hamiltonian, neglecting some terms of Hi 
and i?2 as described above. With the parameters de- 
rived from optical experiments [Eq. (jl6|) ]. we obtain 
t^ff — 0.101 eV, which is in very good agreement with 
the estimate of Rcf. 72. 

The final form of the non-interacting part of the effec- 
tive Hamiltonian contains the trigonal distortion D 0.1 
eV and a smaller direct hopping between Co ions t'.— 
These terms are beyond our calculations based on the 
diagonalization of a CoOe cluster. The former because 
we assumed Oh symmetry and the latter because direct 
Co-Co hopping cannot be included in the cluster. Nev- 
ertheless, the largest energies in the problem are the in- 
teractions given by Eqs. (jl7l) within our approach. Note 



that the value of the trigonal splitting 3-D is an order of 
magnitude smaller than the effective cubic crystal field 
(near 3 eV including the effects of Co-0 hybridization) 
and therefore does not sensibly affect the derivation of 
the remaining parameters of -ffcff, which were obtained 
assuming Oh symmetry. 

E. The effect of the overlap 

While in Hcs one assumes a basis with orthogonal 
states, the states of H involved in the calculation of ^eff 
above are non-orthogonal, due to the fact that some p or- 
bitals belong to different clusters. In the cuprateSfi^ and 
in double perovskites,^^ the problem of non-orthogonality 
has been solved by a change of basis to orthogonal Wan- 
nier functions centered at transition-metal atoms. In the 
present case, in which five different 2p O Wannier func- 
tions centered at each Co site should be considered, this 
procedure cannot be handled analytically and becomes 
very awkward. Fortunately, the larger number of or- 
bitals involved has also the consequence that the overlap 
between many-body states is considerably reduced with 
respect to those of the one-particle Wannier functions. 
This also takes place to a lesser extent in the cuprates, 
in which the overlap between Zhang-Rice states is S=- 
1/8, and leads to a rapidly convergent expansion of Hcs 
in terms of S 

Furthermore, a close examination of the structure of 
Naa;Co02 shows that in all non-zero elements of the over- 
lap matrix between bonding combinations of O 2p or- 
bitals [see Eq. p^ ]. at least one O t2g orbital is involved. 
Since the occupation of the latter is small in the relevant 
many-body eigenstates (see Table I) , the overlap between 
the latter is reduced. In particular we find 

5= (g,(4)5,(5,zx)|g,(5,2;z)5,(4)) - (1|2) = -0.0357, 

(19) 

where for simplicity, we denoted as |1) (|2)), the many- 
body state in which all clusters except j (i) are in the 
ground state for four holes and the remaining cluster is 
in the lowest lying state for 5 holes and symmetry zx 
(jjz). To first order in the overlaps {l\m) — Sim, m ^ I, 
orthonormal states can be obtained using 

10 = 10-^ E ^"'i"^)- (20) 

Then 

Us = (i|ff|2) - Hi2 - \Si2{Hii + F22), (21) 

where Him = {l\H\m). Since by symmetry Hn = H22, 
the above result can be written as 

i,ff = (l|(ff-i/22)|2). (22) 
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Except for the sign, the second member coincides with 
the second member of Eq. (fT5)) . Therefore the resuh 
previously obtained is not modified by terms hnear in 
the overlap. 

IV. SUMMARY AND DISCUSSION 

We have constructed the cell Hamiltonian H), that 
consists of a transition-metal atom M and its neigh- 
boring O atoms. This Hamiltonian is the basis to 
construct a low-energy effective model Hcs in which 
the O atoms are eliminated, using the cell-perturbation 
methodfi^ii^"— ii^iii The Hamiltonian has the same form 
if O is replaced by other elements with p states near the 
Fermi energy, like P or As in the Fe pnictides. An es- 
sential part of the cell Hamiltonian and the full starting 
Hamiltonian H is the interaction between electrons inside 
the d shell Hi when all 10 spin-orbitals are important. 
We have constructed Hj in a basis of orbitals appropriate 
for a cubic, hexagonal, tetragonal or tetrahedral environ- 
ment of the M atoms. For other symmetries a change of 
basis may be required. 

While Hj has a trivial form U ni^riii in the Hub- 
bard model where only one orbital per site i is relevant, 
the interaction contains interorbital Coulomb repulsion, 
exchange and pair hopping terms when only either eg 
[Eq.©] or t2g [Eq.®] are involved. These already makes 
the physics of transition-metal compounds richer than 
that of the Hubbard model^i^.23i37-39^^ 

If all orbitals should be retained, more complicated 
terms appear [Eqs. © and pH]) ]. Our results suggest 
that because of the large M-0 hopping for Cg holes in 
cubic or tetragonal symmetry, they should be retained in 
general, even if the formal configuration of M contains 6 
d electrons or less (corresponding to only t2g electrons in 
a ionic picture). 

In the case of the cobaltates NaxCo02, where mainly 
3d'' and 3d^ configurations play a role, our results as 
well as previous ones,''" show that addition or removal 
of electrons in the d shell causes a strong charge re- 
distribution between metal Cg states and linear combi- 
nation of 2p O states with the same symmetry. This 
has important consequences for the parameters of i?cff- 
In particular, the effective Coulomb repulsion UcS is 
smaller than previously assumed in calculations of HcS 
which include the effect of correlations As a conse- 
quence if the hopping parameters and the trigonal dis- 



tortion D are obtained from a fit of the bands obtained 
from first principles, six hole pockets appear along the 
T — K directions,—"^ which are not detected in pho- 
tocmission experiments."^'^ Instead, using the values of 
the interactions and effective hopping obtained as de- 
scribed in the previous Section, and the larger value of 
D given by quantum-chemistry configuration-interaction 
calculations,^^ these pockets are absent according to cal- 
culations which apply the dynamical mean-field theory 
(DMFT) to Hgffi^ The resulting electronic dispersion 
near the Fermi energy agrees with experiment. Recent 
photoemission experiments in misfit cobaltates show re- 
sults similar as previous ones, with a significant band 
reduction due to correlations 

In general, the local-density approximation (LDA) un- 
derestimates gaps and one-particle excitations energies. 
Thus one might suspect that it underestimates the trig- 
onal distortion energy D in the cobaltates. The above 
results suggest that taking the one-body parameters of 
the effective model HcS with metal sites only from a 
fit of the bands obtained in LDA is not valid in gen- 
eral, or at least when the degree of covalency is impor- 
tant. This is also the case of NiO, for which agree- 
ment with experiment in LDA-I-DMFT calculations is 
only achieved once the O bands are explicitly included 
in the model, or when the O atoms have been inte- 
grated out using low-energy reduction procedures sim- 
ilar as ours, which take into account correlations from 
the beginning^i2i. The research in the superconduct- 
ing cuprates also supports the above statement. In this 
case, the parameters of the effective one-band Hubbard 
OT t — J models are obtained accurately using system- 
atic low-energy reduction procedures^^^— from a multi- 
orbital Cu-0 model^^'^^ with parameters obtained from 
constrained-density- functional theory4^i^ 

For the cobaltates, we have determined the parameters 
of the original multiorbital model form a fit of polarized 
x-ray absorption spectra of Naa:Co02f^ 
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